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The one-loop renormalization in field theories can be formulated in terms of the heat kernel ex- 
pansion. In this paper we calculate leading contributions of discontinuities of background fields and 
their derivatives to the heat kernel coefficients. These results are then used to estimate contribu- 
tions of the discontinuities to the Casimir energy. Sign of such contribution is defined solely by 
the order of discontinuous derivative. We also discuss renormalization in the presence of singular 
(delta-function) potentials. We show that an independent surface tension counterterm is necessary. 
This observation seems to resolve some contradictions in previous calculations. 
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I. INTRODUCTION 

Recent years have seen much progress in understanding the Casimir effect Q, ■ Although this effect is oftenly 
considered as a macroscopic manifestation of quantum fields, many ingredients of quantum field theory have not been 
yet reformulated to the Casimir framework. Renormalization is probably the most important example of such an 
ingredient. Early calculations of the Casimir force between rigid bodies were not affected by the renormalization 
problems since the divergent part of the Casimir energy is local (does not depend on the distance) and, therefore, 
does not contribute to the force. There are, however, many quantities of interest (like the Casimir stress) which do 
contain a divergent part and require some renormalization. 

There exists a subtraction prescription [2|, |3( which allows to obtain well-defined results for the Casimir energy even 
if it is divergent. This prescription consists in expanding the regularised Casimir energy in log-power series for large 
mass m of the fluctuating field and subtracting all non-negative powers of m and logm terms. A similar prescription 
is being used in quantum field theory in curved space Q. In two dimensional scalar theories it is equivalent to the 
"no tadpole" condition p|. Unfortunately, no such statement is known for other models or other dimensions. Besides, 
this large mass prescription is not applicable for massless fields. 

In order to remove ambiguities in the Casimir energy, it was suggested just a couple of years ago 0,0 to perform 
the calculations in renormalizable theories. This direction attracted a lot of interest, but the results of different groups 
are still contradicting (cf. 8, 9] and references therein). 

The main difficulty in applying the renormalization procedure of quantum field theory to the Casimir-type problems 
is that the latter one usually assume singular backgrounds (which may be penetrable walls, or boundary conditions, 
or just non-smooth potentials). This motivated us to study quantum field theory (in the one-loop approximation) in 
the presence of non-smooth background fields. Namely, wc consider the case when pth derivative of the background 
field jumps on a surface E of co-dimension one. 

In this paper we use the zeta function regularization and the heat kernel methods (both are sketched in sec.[n]). In 
this regularization, the one- loop divergences are expressed in terms of the heat kernel coefficients dfc. These coefficients 
are analysed in sec. lllll for non-smooth potentials. In sec. HII Al we calculate the linear order in the potential and find 
how the localised heat kernel is modified by the presence of singularities. Somewhat surprisingly, the global heat kernel 
(which defines the counterterms) is not sensitive to the singularities in this order. The quadratic order in the potential 
is studied in sec. HII Bl where we calculate leading contributions to global heat kernel. There we prove a "folk theorem" 
that if standard "smooth" expression for is divergent due to the singularity, a non-zero surface contribution to 
dk-i should appear. These results allow us to calculate the leading contribution of the singularity to the Casimir 
energy (sec. IIII C|l . Then, in sec. II VI we discuss renormalization of theories with delta-function singularities. 
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II. THE REGULARIZATION 



In this section we sketch basic technical results on the zeta-function regularization and the heat kernel. For a more 
detailed introduction one may consult the monographs or the review ^^l- The reader may also consult a recent 
paper by Fulling ^| which deals with applications of the heat kernel technique to the Casimir energy calculations. 

To simplify our discussion in this paper we consider scalar field theories only. The classical action reads 

c= / M dn *Q(v^) 2 +tW) • a) 

We suppose that the manifold M is fiat, so that it is either torus or W 1 . For M = W l one has to assume some 
fall-off conditions on the fields to obtain convergent integrals. We assume also that the metric on M has Euclidean 
signature. We use the background held formalism, so that we split (p = $ + <j) where $ is a background field, and 
is a quantum fluctuation. Note, that the background field $ may also include a part which describes "external 
conditions" (like, e.g., boundary conditions or domain walls or membranes). To calculate the one- loop effective action 
one should expand the action about the background field $ keeping the quadratic order terms in fluctuations only. 



C 2 = \ f d?xct>D[$\<t>, (2) 

where D[Q] is an operator of Laplace type, 

D = -(\7 2 + E), E = -U"(&). (3) 

For a multicomponent <p the potential E is matrix-valued. 

Formal path integration over <f> leads to the following result for the one-loop effective action 

W = -lndetOD). (4) 

Right hand side of @ is divergent and has to be regularised. To this end we use an integral representation for the 
determinant 



1 f°° dt 



(5) 



with a regularization parameter s (the regularization is removed in the limit s — > 0). The heat kernel K(t,D) is 
defined as a functional trace, 

K{t,D)=Tr L ,{e- tD ). (6) 
Let us define the zeta function of D by the equation: 

(( S ,D) = Tr L2 (D- s ). (7) 
Now we can rewrite W s in terms of the zeta function 

w s = -^ 2s r( s )<;(s,D). (8) 

This procedure is called the zeta-function regularization |10|. 
The regularised effective action JHJ has a pole at s — 0: 

W s = -\(-- lE + \nA C(0, D) - V(0, D) , (9) 



2 \s 



where is the Euler constant. 

There is an asymptotic series as t — > +0 



K(f,t,D) = Tr L2 (/cxp(-t J D)) = £ t^/ 2 a k (f, D) . (10) 



fc>0 
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(here we have introduced a heat kernel smeared with a smooth function /. It is related to the unsmeared kernel 
© by K(t,D) — K(l,t, D).) If E is smooth and M has no boundary, all odd numbered coefficients in IjlUI) vanish, 

«2fe+l = 0. 

There is an important relation 

((0,D) = a n (l,D) (11) 

which tells that the one-loop divergences are defined by the heat kernel asymptotics. We remind that n is dimension of 
the underlying manifold M. Besides that, the heat kernel expansion defines short-distance behaviour of the propagator 
and the large mass expansion of the effective action. 

We shall use also a bi-local kernel K(x,y;t) of e~ tD . It is related to the smeared heat kernel by means of the 
equation 

K(f,D,t)= I dxtrf(x)K(x,x;t). (12) 

J M 

Here tr denotes a matrix trace over all discrete indices if they are in the model. It should be distinguished from Tr^2 
which is a functional trace in the space of square integrable functions. 

Let us consider a surface S of co-dimension one in M. We take a background field $ which is smooth everywhere 
except for S where its' pth normal derivative has a discontinuity. The potential E shall then possess the same property. 
At least locally, S divides M into two parts M + and M~ . Let v + and v~ be unit normal vectors to £ pointing inside 
M + and M~ respectively. Since M is smooth, v~ = but this may be not true for more general geometries. Let 

EK fe )+ (respectively E^~) be a limit of fcth derivative w.r.t. v + (resp. v~) of E(x) as x — > £ from the M + (resp. 
M~) side. If the fcth derivative is continuous, £j( fe )+ = (— l) k E^~ . For example, if M = R 1 and £ = {x = 0}, then 
M ± =R±, E^ = ±d x E. 



III. HEAT KERNEL FOR NON-SMOOTH POTENTIALS 



In this section we restrict ourselves to the case of flat £. For a non-flat £ the heat kernel expansion should contain 
additional terms with extrinsic curvature of £. Such terms have larger canonical (mass) dimensions than the terms 
considered below and shall, therefore, contribute to higher heat kernel coefficients. 



Local heat kernel in the linear order 



To analyse the heat kernel asymptotics we use the perturbative expansion [l6l | (see also |12| for a short overview). 
The exponent exp(— tD) = exp(t(A + E)) with A = can be expanded in a power series in E: 

t t S 2 



e- tD = e tA 



This expansion is purely algebraic. Each order of E in (|13|) is given by a convergent integral if E is smooth or has a 
singularity located on a surface of co-dimension one (this can be a discontinuity of a derivative or even a delta- function 
singularity). In the case of ^-singularities on a submanifold of co-dimension two or higher, there might be problems 
with the convergence. More careful estimates can be found in Ref. |l7j . 

Here we analyse the heat kernel expansion to the linear order in E. Let / be a smooth function on M. Let 
Ko(x,y,t) be the heat kernel for E = 0. Then to this order 

K(f, E, t) = K(f, 0, t) - [ dr [ dx [ dyf(x)K (x, y,r- t)E(y)K (y, x,t) + ... (14) 
Jo Jm Jm 

Obviously, the second term in 114|) is symmetric w.r.t. exchanging the role of E and /. Therefore, it can be interpreted 
as a linear order term of the heat kernel with the potential / and the smearing function E. This term then reads: 

dxK 1 (x,x,t)E(x) (15) 

M 
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where K\ is the heat kernel with the potential f(x) at the linear order in /. The crucial point is that / is smooth. If 
we neglect all curvatures, then local heat kernel coefficients corresponding to K are: 

a 2 k{x) = (4^)-"/ 2 a 2 (fc)A' £ - 1 /^) (16) 

where [l4^ : 

2fe' 

a >™ = J2ky. < 17 > 
The heat kernel coefficients we are interested in read: 

a 2 fe = (4^)-"/ 2 a 2 (fc) f dxE{x)A k - 1 f(x) (18) 

J M 

Odd- numbered coefficients are all zero. Formula Q18|) can be checked for k = 1, 2 [l5| . 

An interesting observation regarding (|18|l is that all these coefficients vanish for / = 1, i.e. one loop divergences 
and the large mass expansion in the effective action are not affected by non-smoothness of the potential (to the linear 
order studied in this section). 



B. Global heat kernel to quadratic order 

In this section we study the heat trace asymptotics in the order E 2 and prove for this case a statement which existed 
in the folklore for many years: If due to a singularity of the background the volume contribution to a 2 k diverges, the 
coefficient a^k-x should have a non-zero surface contribution. This statement is based on particular case calculations 
for scalar backgrounds 0, 0] and for singular magnetic fields J3 , and on analytic calculations of 0, 0] . A similar 
conclusion for dielectric problems follows from the analysis of [3j . In the present paper we give explicit expressions for 
the leading contribution to the heat kernel expansion from discontinuities of derivatives of the potential of arbitrary 
order. 

We start with analysing relevant volume and surface invariants. For a smooth potential E there is an invariant 
(VP +1 E) 2 (which means (VAp/ 2 £) 2 for p even and (A^+^^E) 2 for p odd). It has dimension (2p + 6) and can 
contribute to the coefficient a 2p +6: 



(47r)-"/ 2 /3(p)tr f dx{V p+1 E) 2 . (19) 

J M 

If pth derivative of E is discontinuous the expression (|19f) contains a delta-function squared and is, therefore, mean- 
ingless. Then we expect that the following heat kernel coefficient appears 

(4 7r )-('»-i)/2 7 (p)tr / dx(E^+ - {-IfE^-f . (20) 



&2p+b 

No other invariant of the same dimension can appear. For example, - {-1)PE^-)(E^+ + (-1)peW~) 

changes sign if one exchanges the role M~ and M~, (E^ + + (-l)PE^-) 2 would give a non-zero 02^+5 even for 
smooth potentials. According to the general theory both f3(p) and 7(p) do not depend on n. We have to make sure 
that they are non-zero. 

We use again the perturbative expansion l|13fl . The heat trace can be also expanded, 

oo 

K(t,D)=Tr(e- tD )=J2KM (21) 
where Kj contains the jth power of E. Localised version of K\ has been studied above. In the next order we have: 

t S2 

ds 2 J ds 1 K (z 7 y;t - s 2 + s 1 )E(y) 
o 

xK (y,z;s 2 -s 1 )E(z). (22) 
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To derive 12211 we used cyclic property of the functional trace in order to combine e ( t ~ s ^ A with e SlA . 

Since the coefficients (3{p) and ^{p) are universal constants we can use some particularly simple model to calculate 
their values. We take M = R, E(x < 0) = 0, E(x) smooth for x > and decreases exponentially fast as x — > +oo. 
We also suppose that first p — 1 derivatives are continuous at x = 0. Then 



K (x,y;t) = (47rfr 1/2 exp 



(x - yf 

At 



We adopt the strategy of |2jj. After removing a redundant integration (|22[1 takes the form: 



K 2 {t) 



t 



^ Jo 

Next we integrate over r to obtain 



dx 



ft exp 
dy I dr 



K 2 (x,y : t) = — / dx dyerfc 



\x - y\ 

Vt 



E(x)E(y) 



E(x)E(y) . 



Let us change the variables 



x = z + rvc, 
2/ = z + r\/t. 



y = z for x > y 
x = z for x < y. 



In both cases fe,r € [0, +oo[. Then 

^3/2 /-oo (-oo 

K 2 {x,y;t) = — dz dr erfc(r)E(z)E(z + rVi) 
4 Jo Jo 

The t — ► asymptotic expansion is now performed by using the following formula: 

/""/(rVi)erfc(r)~fv 



i/2 r ( 1 + ?) f faV n x 

(n+l)!V^ J W 



Now we only have to pick up relevant terms in the expansion. The term with n = 2(p + 1) reads: 

i p+3 1 2fc! 



(47rf)!/2 2 (2fc)! y 
which is consistent with (|19|l for smooth potentials if 



dzE{z)d 2 z {v+1) E{z) , 



p(p) = (-ir +1 -a 2 (p + 2). 

a 2 is given by (|17|l above. Equation l|3U|) follows easily from [l4| and can be used as a consistency check. 
The term with n = 2p + 1 has the form: 

/>oo 

t P+2 4 -( P +2) / dzE(z)dlP +l E{z). 



If first p — 1 derivatives of E vanish at z = 0, this result is consistent with (|2(J|I and gives 

lp = (_i)p+i 2 - 2 p- 5 . 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



We see, that both constants 7 and @ are non-zero. 

Our calculation also confirms that as long as general "smooth" formulae for the heat kernel coefficients give conver- 
gent integrals no modfications appear due to the singularities. In the present context this means that old non-modified 
formulae are valid for ak with k < 2p + 5. 
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C. Application: Casimir energy density at the singularity 

Let us calculate leading contribution of non-smoothness of the potential to the vacuum energy. In the zeta function 
regularization the ground state energy is defined as 

£( s ) = ^ 2fl £4~ 2fl > (33) 

k 

where the regularization parameter s should be taken zero after the calculations, are eigenfrequencies of elementary 
excitations defined as square root of eigenvalues of the Hamiltonian: 

H = -df + V(x) + m 2 = D + m 2 . (34) 

We shall work in a 3 + 1 dimensional theory. Therefore, the operator D is three dimensional. V is a static potential 
which which is supposed to be non-negative. We can rewrite (|33H though the zeta function of H, 



£(s) = ^ 2s Tr (if*-) = l^C (s \,H\ (35) 



The zeta function, in turn, can be expressed through the heat kernel: 



d > - \,H) = r(a 1 _i ) J™ dtt°- 3 / 2 K(t,H) (36) 



and 

K(t, H) = K(t, D)e~ tm2 . (37) 

It is easy to see that the large mass expansion of the vacuum energy is generated by the small t asymptotics of 
K(t, D). Indeed, by substituting the heat kernel expansion in (|35|) with (|36|) and integrating over t we obtain: 

£{s) = l» 2s T ( /_ J2a k (D)m- 2 °- k+i r (« + | - 2) • (38) 

As expected, the terms with k = 0, . . . , 4 are divergent in the limit s — > 0. 

We would like to separate contributions from non-smooth parts of the potential. If pth derivative jumps, the 
leading contribution comes from the coefficient 02^+5 in the heat kernel expansion. Note, that in four dimensions 
all contributions of these type are not divergent, so that we can put s — already in l|38|) . Then the term we are 
interested in reads 

£ lp] = -^«2 P+ 5m- 2rl r (p + ~) • (39) 

Next we use (O and (J22JI to obtain 

£ \p] = (_ 1 )p ( 2 ^il!: 2-3P-7 m -2 P -i I dx ( 6v ( P )y f (40) 

where SV^ is discontinuity of of pth normal derivative of V on S, according to our conventions (—1)!! := 1. 

We see, that discontinuities in the potential itself and in its' even order derivatives tend to increase the vacuum 
energy, while discontinuities in odd order derivatives tend to decrease the vacuum energy. Of course, practically it 
may be not easy to separate contributions from continuous and discontinuous parts. 



IV. RENORMALIZATION WITH SINGULARITIES 



As we have already seen in the presence of singularities the heat kernel expansion is modified. This means that new 
counterterms may appear in quantum field theory at one loop. Non-smoothness of the potential modifies global heat 
kernel coefficients starting with 05 (discontinuous potentials) or even higher (discontinuous derivatives). Therefore, 
they have no effect on counterterms in four dimensions. In order to be closer to physical applications we consider a 
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stronger singularity (delta-potentials) in four dimensions. This example is of particular interest because of extensive 
discussion in the literature (cf. recent works [3, l| and references therein) . 

We start with some technical information. Consider an operator which has a singular part concentrated on the 
surface E: 



D 



-(d 2 +E(x)+v(x)6v(x)). 



(41) 



For simplicity, we do not consider here any gauge fields and suppose that E{x) is smooth. A mathematically correct 
formulation of spectral problem for operators with delta-like singularities yields a spectral problem for the regular 
part of D outside of E supplemented by matching conditions on E. 



IE = 
(1) + 



,d)- 



o. 



(42) 



We remind, that according to our notations is a feth normal derivative of <j> calculated on M + or M side of E. 

For a smooth (f>: 4>^ + = —<j)^~ (cf. sec llTj) . 

The heat kernel coefficient which is responsible for one- loop divergences in four dimensions reads |l7) : 



a 4 (£>) = (47r)-"/ 2 tr 



dx-E 2 + 



A I 



Ev 



(43) 



Note, that no additional terms appear in (|43[1 even if the surface E is curved. This implies that one needs the same 
counterterms on spherical surfaces as on flat ones. 

Consider now an action containing a surface interaction term: 



C 



(44) 



Here m, A, A, er, and A are some coupling constants. Next we introduce a continuous 1 background field ip = $ + <f> 
and expand the action 1)44(1 up to quadratic order in quantum fluctuations <fi. As a result, we obtain (j2J where D is 
given by (|4*T|l and 



E = -m z - A$ 2 



v = —X. 



(45) 



With the help of (JSJ, l(TT)l. l(4^|) and lUSTl we immediately obtain that the one loop divergences can be cancelled by 
the following counterterms: 



d .\ = l^L +0 ( S °l Sm 2 =l^- + O(s ), 5 X=-^ + 0( S °), 



s 32tt 2 



S 167T 2 



S 167T 2 



firr = IJ— (-H 3 + Xm 2 ) + O(s ), SX = + O(s ), 



s 32tt 2 V 6 



s 32tt 2 



(46) 



We see that the model (|44f) is renormalizable at least at one loop. 

It is important that both constants A (the cosmological constant) and a (the surface tension) must be present in 
the action to achieve renormalizability. One can consistently put A = to all orders of perturbation theory since its' 
observed value is negligible. The surface tension a cannot be excluded on the same grounds. It has to be considered 
as an (experimental) input, as other coupling constants. To remove all ambiguities in the model one needs five 
independent normalisation conditions (including the "trivial" one A = 0) which should fix the O(s ) terms in l|4t)|l . 
To fix Sm 2 and 5A one can proceed as in non-singular theories. Namely, one can consider the effective action for a 
slowly varying background field $ located far away from E. Then Sm 2 and A can be determined by prescribing certain 
values to the second and fourth derivatives of the effective action w.r.t. <& at $ = 0. Physically this is equivalent 
to fixing position of the pole in the propagator of bosons and to fixing the value of the four-boson vertex for zero 
external momenta (both processes have to be considered very far from E). In principle, finite renormalization of 
A can be fixed by relating the renormalized value of A to some amplitudes of scattering of ^ on H. It is not clear 
however which condition is the most convenient one. Unless all these conditions have been formulated one cannot give 



1 A stronger requirement that <J> is smooth would exclude classical solution of this model. 
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a scheme-independent meaning to the surface tension. This provides an alternative explanation to the contradictions 
reported in We should also note, that if one is interested only in the one-loop Casimir energy and stress on the 
background with $ = (as in @), ambiguities in finite renormalization of to 2 , A and A are not important since they 
are of the order h and do not enter the classical energy for $ = 0. 

It is instructive to compare (|44|l to other models which appeared in the literature. The model considered by Milton 
did not contain the self-interaction A (which is not important) and the surface tension a (which is very important 
for discussing the counterterms) . The model of the MIT group Q uses coupling to an external source. Although, 
the counterterm action suggested in [j| indeed allows to remove all divergences (at least as long as the source is 
smooth), the model shares a common difficulty of all models with non-dynamical external fields: the latter either have 
to be fixed artificially or require infinite number of "experiments" to be properly determined (if one "experiment" 
allows to determine one real number). Indeed, each coefficient in the Taylor expansion of the external source (or each 
form-factor) plays the role of an independent coupling. 

There is a somewhat exotic example of a model with a surface interaction, 

C = / A /*Q( V ^) 2 + ^ 4 ) + j/ x {v\ (47) 

which does not require any surface tension counterterms in the zeta-regularization at one loop. In other schemes 
(momentum cut-off, e.g.) such counterterms may appear also for (|47|l . 



V. CONCLUSIONS 



In this paper we studied quantum field theory on non-smooth backgrounds. We have analysed leading contributions 
of the singularities to the heat kernel coefficients. We found that influence of the non-smoothness is rather mild: no 
new counterterms appear in dimension up to 4. We were also able to prove an important statement on the behaviour 
of the heat kernel expansion for singular backgrounds: "smooth" formulae for the heat kernel coefficients are valid as 
long as they are convergent. If the "smooth" expression for aik diverges for some fc, a non-zero surface term in a2k-i 
inevitably appears. We demonstrated that discontinuities of even order derivatives contribute a positive amount to 
the Casimir energy, while contributions from discontinuities of odd order derivatives are negative. 

We have also analysed renormalization of theories with stronger (delta-function) singularities. Our main message 
is that an independent counterterm to the surface tension is needed in this case. Therefore, the one loop surface 
tension crucially depends on the normalisation condition. If no such condition is specified, one cannot give a scheme- 
independent meaning to the surface tension. 

One may extend our results by using the heat kernel expansion with non-smooth gauge fields |2lj |. non-smooth 
geometries 0, other surface singularities [2i|, and, perhaps, even singularities of conical type p3[ . 
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